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As a basis for the interpretation of data gathered by LOPES and other experiments, 
we have carried out Monte Carlo simulations of geosynchrotron radio emission 
from cosmic ray air showers. The simulations, having been verified carefully with 
analytical calculations, reveal a wealth of information on the characteristics of 
the radio signal and their dependence on specific air shower parameters. In this 
article, we review the spatial characteristics of the radio emission, its predicted 
frequency spectrum and its dependence on important air shower parameters such 
as the shower zenith angle, the primary particle energy and the depth of the shower 
maximum, which can in turn be related to the nature of the primary particle. 



1. The simulations 

Two main mechanisms have been considered to contribute to radio emis- 
sion from cosmic ray air showers: Askaryan-type Cerenkov radiation arising 
from a negative charge excess moving through the atmosphere at velocities 
faster than the speed of light in air, and emission generated as a conse- 
quence of the deflection of charged particles in the earth's magnetic field. 
A number of historical results illustrated that, while the Cerenkov emis- 
sion mechanism dominates in dense media such as ice, the geomagnetic 
mechanism is dominant in the atmosphere. Our simulations thus focus on 
the latter mechanism, interpreting the radio emission as "coherent geosyn- 
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chrotron radiation" arising from the geomagnetic deflection of highly rela- 
tivistic electron-positron pairs generated in the air shower cascade 1 . 

In order to gain a solid understanding of the emission characteristics, in 
a first step, we have performed analytical calculations of the expected radio 
signal 2 . In a second step, we have then improved on the analytical results 
with detailed Monte Carlo simulations, which we have directly compared 
and thus verified with the analytical results and the available historical 
data 3 . While these Monte Carlo simulations still incorporate a somewhat 
simplified air shower model based on analytical parametrizations, they do 
already take into account the most important air shower characteristics 
such as longitudinal (arrival time) and lateral particle distributions, energy 
and track-length distributions, the overall longitudinal development of the 
air shower and the geometry of the air shower and magnetic field. Our 
simulations thus for the first time have provided a prediction of the radio 
emission from realistically modeled air showers. Furthermore, the model is 
currently being enhanced by substituting the parametrized particle distri- 
butions with CORSIKA 4 -generated distributions. 

We here present only a subset of the results derived with our Monte 
Carlo code. A more detailed analysis, including a parametrization of the 
derived dependences, has been published elsewhere 5 . 

2. Simulation results 

We first consider the very simple scenario of a vertical air shower with pri- 
mary particle energy of 10 17 eV developing to its maximum at ~ 630 g cm~ 2 , 
i.e., at ~ 4000m above ground. The geomagnetic field is adopted as 70° 
inclined with a strength of 0.5 G, corresponding approximately to the field 
configuration in central Europe. 

Figure 1 shows the ground-level total field strength emission pattern at 
10 MHz, visualized as a contour plot, and the frequency spectra derived at 
various radial distances from the shower center. The emission pattern shows 
remarkable symmetry and is almost circular. This is not a trivial result, as 
the emission process itself, i.e., the deflection of electrons and positrons in 
the geomagnetic field, is a highly directed process. The circularity of the 
footprint illustrates that most of the emission stems from particles having 
short track lengths. A slight north-south asymmetry introduced by the 
inclination of the geomagnetic field is also visible. The frequency spectra 
shown in the right panel illustrate that the field strength drops quickly to 
higher observing frequencies. This is a direct consequence of diminishing 
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Figure 1. Radio emission from a vertical 10 eV air shower. Left: 10 MHz total field 
strength emission pattern. Right: Frequency spectra at (from top to bottom) 20 m, 
140 m, 260 m, 380 m and 500 m north of the shower center. 



coherence as the wavelength becomes shorter and thus comparable to the 
dimensions of the air shower pancake, in particular its thickness of a few 
meters. The decrease is stronger at larger distances from the shower center. 
When one enters the incoherent regime the frequency spectra exhibit un- 
physical seeming features such as rapidly alternating series of maxima and 
minima. Realistic calculations of the emission in this regime can only be 
performed with a better underlying air shower model taking into account 
inhomogeneities in the shower in very mildly thinned calculations. 

Figure 2 demonstrates the changes arising in the transition from a verti- 
cal to a 45° inclined air shower (coming from the south). The emission pat- 
tern becomes elongated considerably along the shower axis. This is mainly 
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Figure 2. Radio emission from a 45° inclined 10 17 eV air shower. Left: 10 MHz total 
field strength emission pattern. Right: Frequency spectra at (from top to bottom) 20 m, 
140 m, 260 m, 380 m and 500 m north of the shower center. 
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a projection effect directly associated with the inclination of the shower 
axis. On closer look, however, the emission pattern becomes wider (and 
less peaked) as a whole, even in the direction perpendicular to the shower 
axis. The reason for this is that the maximum of the inclined shower lo- 
cated at the same (slant) atmospheric depth of ~ 630gcm~ 2 now is at 
much greater geometrical distance from the observer at ground-level. This 
geometric effect has direct influence on the slope of the radio emission's lat- 
eral distribution. A look at the frequency spectra in the right panel shows 
that coherence is also retained up to higher frequencies in case of inclined 
showers. Their larger radio footprint combined with the large solid angle 
associated with medium to high zenith angles thus makes inclined showers a 
particularly interesting target for radio observations 5 ' 6 . At near horizontal 
inclination, even neutrino-induced air showers might become observable. 

Figure 3 illustrates two additional parameters that have direct influ- 
ence on the radio signal. The left panel shows the impact of the primary 
particle's energy. The electric held strength at all distances scales as a 
power-law with the primary particle energy. The power-law index is very 
close to unity, i.e., that of the linear relation expected for coherent emis- 
sion. To larger distances, the slope of the power-law gets flatter due to 
the effect that more energetic showers on average penetrate deeper into the 
atmosphere and thus have their shower maximum geometrically closer to 
the observer. As already discussed in the context of inclined showers and 
illustrated in the right panel, this directly influences the lateral distribution 
of the radio emission. Since the depth of the shower maximum can in turn 
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Figure 3. Radio emission from vertical air showers. Left: Scaling of the 10 MHz electric 
field strength with primary particle energy at various distances from the shower center. 
Right: Dependence of the radio signal's lateral distribution on the depth of shower 
maximum X max - -fmax in g cm -2 is for red/solid: 560, green/dashed: 595, blue/dotted: 
631, violet/short dotted: 665, turquois/dash-dotted: 700, black/double-dotted: 735. 
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be related to the nature of the primary particle, its influence on the radio 
signal's lateral distribution can potentially be used to probe the primary 
particle composition with radio measurements 7 . 

Another important result of the simulations (not shown here explicitly) 
are the predicted linear polarization characteristics of the radio signal 5 . 
They can be used to directly verify the geomagnetic origin of the emis- 
sion. To make the simulation results available for easy comparison with 
experimental data, they are also available as a parametrization formula 5 . 

3. Conclusions 

We have carried out elaborate Monte Carlo simulations of geosynchrotron 
radio emission from cosmic ray air showers. Special care has been taken to 
verify the Monte Carlo results with analytical calculations and historical 
data, giving us a good understanding of the emission process and thus solid 
confidence in the predictions. The simulations predict many important 
characteristics of the radio emission and their relation to parameters of the 
associated air shower. The total field strength emission pattern is very 
regular and symmetric in the coherent regime. The geomagnetic origin of 
the emission can be directly verified with polarization measurements. The 
frequency spectra cut off quickly to high frequencies, making low observing 
frequencies around a few tens of MHz desirable. Inclined air showers have 
a much wider emission pattern and are thus particularly suitable for radio 
observations. The slope of the lateral distribution can be directly related 
to the geometrical distance between observer and shower maximum. It is 
thus not only sensitive to the shower zenith angle, but also to the nature of 
the primary particle. As expected for coherent emission, the electric field 
strength scales approximately linearly with the primary particle energy. 
These predictions will allow to analyze and interpret experimental data 
such as those provided by LOPES 8 and other experiments. 
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